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Abstract 

Consideration is given to thermodynamical properties of a three-dimensional Bose- 
condensate of translation-invariant bipolarons (Tl-bipolarons). The critical temperature 
of transition, energy, heat capacity and the transition heat of Tl-bipolaron gas are calcu¬ 
lated. The results obtained are used to explain experiments on high-temperature super¬ 
conductors. 
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1 Introduction. 

The theory of superconductivity is one of the hnest and old subject-maters of condensed mat¬ 
ter physics which involves both macroscopic and microscopic theories as well as derivation of 
macroscopic equations of the theory from microscopic description [1] . In this sense the theory 
was thought to be basically completed and its further development was to have been concerned 
with further detalization and consideration of various special cases. 

The situation changed when the high-temperature superconductivity (HTSC) was discov¬ 
ered [2]. Surprisingly, it was found that in oxide ceramics, the correlation length is some orders 
of magnitude less than that in traditional metal superconductors while the ratio of the energy 
gap to the temperature of superconducting transition is much greater [3]. The current status 
of research can be found in books and reviews a-na. 

Today the main problem in this field is to develop a microscopic theory capable of explaining 
experimental facts which cannot be accounted for by the standard BCS theory. One might 
expect that development of such a theory would not affect the macroscopic theory based on 
phenomenological approach. 
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With all the variety of modern versions of HTSC microscopic descriptions: phonon, plasmon, 
spin, exciton and other mechanisms, the central point of constructing the microscopic theory is 
the effect of electron pairing (Cooper effect). In what follows such a bosonization of electrons 
provides the basis for the description of their superconducting condensate. 

The phenomenon of pairing, in a broad sense, is considered as arising of bielectron states, 
while in a narrow sense, if the description is based on phonon mechanism, it is treated as 
formation of bipolaron states. For a long time this view was hindered by a large correlation 
length or the size of Cooper pairs in BCS theory. For the same reason, over a long period, the 
superconductivity was not viewed as a boson condensate (see footnote at p. 1177 in [E]). A 
signihcant reason of this lack of understanding was a standard idea that bipolarons are very 
compact particles. 

The most dramatic illustration is the use of the small-radius bipolaron (SRB) theory to 
describe HTSC [10], [IT], [18]. It implies that a stable bound bipolaron state is formed at one 
node of the lattice and subsequently such small-radius bipolarons are considered as a gas of 
charged bosons (as a variant individual SRP are formed and then are considered within BCS of 
creation of the bosonic states). Despite the elegance of such a picture, its actual realization for 
HTSC comes up against inextricable difficulties caused by impossibility to meet antagonistic 
requirements. On the one hand, the constant of electron-phonon interaction (EPI) should be 
large for bipolaron states of small radius to form. On the other hand, it should be small for 
the bipolaron mass (on which the superconducting temperature depends [in]-[21]) to be small 
too. Obviously, the HTSC theory based on SRP concept which uses any other (non-phonon) 
interaction mechanisms mentioned above will run into the same problems. 

Alternatively, in describing HTSC one can believe that the role of a fundamental charged 
boson particle can be played by large-radius bipolarons (LRB) [25]-[29]. Historically just this 
assumption was made by Ogg and Schafroth I2SI, [30] long before the development of the 
SRP theory. When viewing Cooper pairs as a peculiar kind of large-radius bipolaron states, 
one might expect that the LRP theory should be used to solve the HTSC problem. 

As pointed out above, the main obstacle to consistent use of the LRP theory for explaining 
high-temperature superconductivity was an idea that electron pairs are localized in a small 
region, the constant of electron-phonon coupling should be large and, as a consequence, the 
effective mass of electron pairs should be large. 
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In the light of the latest advances in the theory of LRP and LRB, namely, in view of 
development of an all-new concept of delocalized polaron and bipolaron states translation- 
invariant polarons (Tl-polarons) and bipolarons (Tl-bipolarons) [31]-[37], it seems appropriate 
to consider their role in the HTSC theory in a new angle. 

We recall the main results of the theory of Tl-polarons and bipolarons obtained in pT]-[37j. 
Notice that consideration of just electron-phonon interaction is not essential for the theory and 
can be generalized to any type of interaction. 

In what follows we will deal only with the main points of the theory important for the 
HTSC theory. The main result of papers [3T] - 1371 is construction of delocalized polaron and 
bipolaron states in the limit of strong electron-phonon interaction. The theory of Tl-bipolarons 
is based on the theory of Tl-polarons [21], [32] and retains the validity of basic statements 
proved for Tl-polarons. The chief of them is the theorem of analytic properties of the ground 
state of a Tl-polaron (accordingly Tl-bipolaron) depending on the constant of electron-phonon 
interaction a. The main implication of this statement is the absence of a critical value of 
the EPI constant Oc, below which the bipolaron state becomes impossible since it decays into 
independent polaron states. In other words, if there exists a value of ac, at which the Tl-state 
becomes energetically disadvantageous with respect to its decay into individual polarons, then 
nothing occurs at this point but for a < ac the state becomes metastable. Hence, over the 
whole range of a variation we can consider Tl-polarons as charged bosons capable of forming 
a superconducting condensate. 

Another important property of Tl-bipolarons is the possibility of changing the correlation 
length over the whole range of [0, cxd] depending on the Hamiltonian parameters [31]. Hence, 
it can be both much larger (as is the case in metals) and much less than the characteristic size 
between the electrons in an electron gas (as happens with ceramics). 

A detailed description of the theory of Tl-polarons and bipolarons and description of their 
various properties is given in review [37] . 

An outstandingly important property of Tl-polarons and bipolarons is the availability of an 
energy gap between their ground and excited states (§3). 

The above-indicated characteristics can be used to develop a microscopic HTSC theory on 
the basis of Tl-bipolarons. 

The paper is arranged as follows. In §2 we take Pekar-Froehlich Hamiltonian for a bipolaron 
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as an initial Hamiltonian. The results of three canonical transformations, such as Heisenberg 
transformation, Lee-Low-Pines transformation and that of Bogolyubov-Tyablikov are briefly 
outlined. Equations determining the Tl-bipolaron spectrum are derived. 

In §3 we analyze solutions of the equations for the Tl-bipolaron spectrum. It is shown 
that the spectrum has a gap separating the ground state of a Tl-bipolaron from its excited 
states which form a quasicontinuous spectrum. The concept of an ideal gas of Tl-bipolarons is 
substantiated. 

With the use of the spectrum obtained, in §4 we consider thermodynamic characteristics of 
an ideal gas of Tl-bipolarons. For various values of the parameters, namely phonon frequencies, 
we calculate the values of critical temperatures of Bose condensation, latent heat of transition 
into the condensed state, heat capacity and heat capacity jumps at the point of transition. 

In §5 we discuss the nature of current states in Bose-condensate of Tl-bipolarons. It is 
shown that the transition from a currentless state to a current one is sharp. 

In §6 the results obtained are compared with the experiment. 

In §7 we consider the problems of expanding the theory which would enable one to make a 
more detailed comparison with experimental data on HTSC materials. 

In §8 we sum up the results obtained. 

2 Pekar-Proehlich Hamiltonian. Canonical transforma¬ 
tions. 

Following [33] - 1371, in describing bipolarons we will proceed from Pekar-Froehlich Hamiltonian: 

^ ^ + U{\n-f2)\+ ( 1 ) 

k 

V + H.c) , U {\n - rll) = .f 

“V J eoolri-ral 

where fq rq are coordinates of the hrst and second electrons, respectively; a^, are operators 

of the birth and annihilation of the held quanta with energy = hwo] m* is the electron 

effective mass; the quantity U describes Coulomb repulsion between the electrons; 14 is the 

function of the wave vector k\ 

e l2'Khvo tvjjQ {1 e^u 

- 4 \ ( 2 ) 
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where e is the electron charge; €oo and Cq are high-frequency and static dielectric permittivities; 
a is the constant of electron-phonon interaction; V is the systems volume. 

In the system of the center of mass Hamiltonian ([1]) takes the form: 


^2 fc2 kr / -- \ 

^ + E + f^(k1) + E 2^fcCOSy + i/.c. j , 

® fe k 

R = (ri + ^ 2 )/‘^i r = fi— f 2 , Me = 2m*, /Xe = m*/2 (3) 

In what follows in this section we will believe h = 1, a;° = 1, Me = 1 (accordingly /ig = 1/4). 

The coordinates of the center of mass R can be excluded from Hamiltonian ([3]) using Heisen¬ 
bergs canonical transformation [38] : 

Si = exp <1 —i kal |> R, S^^akSi = S^^a'^Si = (4) 


Accordingly, the transformed Hamiltonian will be written as: 


Uj^ 1 / _ ^ \ 

H = Si^HSi = -2Ar + U{\f\) + -I- ^ 214 cosy(a/, + a^) + - I ^ka'^at J (5) 

k k \ k J 

From ([3]) it follows that the exact solution of the bipolaron function is determined by the wave 
function '^(r), which contains only relative coordinates r and, therefore, is translation-invariant. 
Averaging of H over 'il>{r) yields the Hamiltonian H: 
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H = 


V fc / k 


hfc(ofc -|- ) -|- T + U, 

k 


% = 2Vk{^ 


kf 


cos- 


v^), f/= (T|f/(r)|T), f = -2(T|A,|T). (6) 


Eq. ([6]) suggests that the bipolaron Hamiltonian differs from the polaron one in that in the 
latter the quantity 14 is replaced by 14 and the constants T, U are added. 

With the use of Lee-Low-Pines canonical transformation [SH] : 


52 



(7) 


where fk are variational parameters having the sense of the distance by which the held oscillators 
are displaced from their equilibrium positions: 


S 2 — flfc T fky S 2 0,1. Sj. — Oj^ fi~, 


( 8 ) 
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for Hamiltonian H: 


(9) 


H = S^^HS2, 

we get: 

H = Ho + i/i, 


Ho = 2j2Vkh + Y,ft + l('E^~f'^ +Wo + f + f/. 

k k \ k / 


T^O — ^ ^ “ 1 “ ^ ^ kk fkfk'i^^k^k' “ 1 “ “ 1 “ ^^k' 


/c,fc' 




1 + —+ fc^fc'/fc'- (10) 


Hamiltonian Hi contains terms linear, threefold and fourfold in the birth and annihilation 
operators. Its explicit form is given in [3T]-[33]. 

Then, as is shown in [ni]-|32], the use of Bogolyubov-Tyablikov canonical transformation 
m for passing on from operators a^, to new operators a^,ak. 


ttk = ^ Mikk’Oik' + ^ ^2kk'(^k' 
k' k' 

^ ^ M2kk'Otk' (11) 

k' k' 

(in which Ho is a diagonal operator), makes mathematical expectation of Hi equal to zero. 

In the new operators a^, ak Hamiltonian fllUl) takes on the form H: 

H = Ehp + 

k 

E,p = AEr + 2Y,%fk + Y,fl + T + U, ( 12 ) 

k k 

where AEr is the so-called recoil energy. The general expression for AEr = AEr {fk} was 
obtained in [32] • Actually, calculation of the ground state energy E^p was performed in [3S| 
by minimization of ffT2|l in fk and in fj. 

Notice that in the theory of a polaron with broken symmetry a diagonalized electron-phonon 
Hamiltonian has the form of ffT2]) |41j . This Hamiltonian can be treated as a Hamiltonian of 
a polaron and a system of its associated renormalized real phonons or as a Hamiltonian whose 
quasiparticle excitations spectrum is determined by flT2|) [42] . In the latter case excited states 
of a polaron are Fermi quasiparticles. 
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In the case of a bipolaron the situation is qualitatively different since a bipolaron is a boson 
quasiparticle whose spectrum is determined by (IT^ . Obviously, a gas of such quasiparticles 
can experience Bose-Einstein condensation (BEG). Treatment of flT^ as a bipolaron and its 
associated renormalized phonons does not prevent their BEG since maintenance of the number 
of particles required in this case takes place automatically due to commutation of the total 
number of renormalized phonons with Hamiltonian flT^ . 

Renormalized frequencies involved in flT^ . according to [31], [32] are determined by the 
equation for s: 



solutions of which yield the spectrum of s = {z/|} values. 


(13) 


3 Energy spectrum of a Tl-bipolaron. 


Hamiltonian flT2]) is conveniently presented in the form: 


n=0,l,2,... 


En — 


Ebpi n 0, 
r'n = Eijp + 0)0 + ^, n 7 ^ 0; 
where in the case of a three-dimensional ionic crystal: 


27r(ni-l) Na- 

Kn = ±-TT-, Ui = 1, -hi, t = 

iV„ ’ * ’ ’ ’ 2 ’ ’ 


(14) 

(15) 


(16) 


Na- - is the number of atoms along the i-th crystallographic axis. 

Let us prove the validity of the expression for the spectrum ([ED, ([ED- Since operators 
an obey Bose commutation relations: 


[ctn, O'na^/ a^/an 

they can be considered to be operators of birth and annihilation of Tl-bipolarons. The energy 
spectrum of Tl-bipolarons, according to f[T3D . is determined by the equation: 


F{s) = 1 


(18) 


where: 



(19) 
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It is convenient to solve equation flT5]l graphically (Fig.l): 



Fig.l. suggests that the frequencies z/fc„ (index i is omitted) lie between the frequencies 
LOk^ and ujkn+i- Hence, the spectrum well as the spectrum are quazicontinuous: 

Ukn — which just proves the validity of flT^ . ffT5|l . 

It follows that the spectrum of a Tl-bipolaron has a gap between the ground state E^p and 
the quasicontinuous spectrum, equal to ojq. 

Below we will consider the case of low concentration of Tl-bipolarons in a crystal. Then 
they can adequately be considered as an ideal Bose gas, whose properties are determined by 
Hamiltonian fll4p . 


4 Statistical thermodynamics of low-density Tl-bipolaron 
gas. 


Let us consider an ideal Bose-gas of Tl-bipolarons which represents a system of N particles 
occurring in some volume V. Let us write No for the number of particles in the lower one- 
particle state and N' for the number of particles in higher states. Then: 


N 


rrin 

n=0,l,2,... 
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1 

g(-En—//)/T — ’ 


n 


( 20 ) 










( 21 ) 


N = N„ + N', N„ = 


g(-Ebp-/i)/T 


n' = j 2 


n^O 


^{En 


In expression N' (HH), we will perform integration over quasicontinuous spectrum (instead of 
summation) ffT^ . ffT5|) and assume /i = E^p. As a result, from ([20]), fl2l]l we get an equation for 
determining the critical temperature of Bose-condensation Tc'. 


Cbp = U (Tc) , (22) 

/. (r.) = {o/t) . F3/.(a) = ^ 

_ / n^l^2TTh^ \ _ _ Wo ^ _ Tc 

\ MeU* ) ’ ^ W* ’ W* ’ 

where n = N/V. The relation between the notation F 3/2 with other notations is given in 
the Appendix. Fig.2 shows a graphical solution of equation fl22]) for the values of parameters 
Me = 2m* = 2mo, where mo is the mass of a free electron in vacuum, w* = 5 meV (~58K), 
n = 10^^ cm“^ and the values: Wi = 0, 2; Cj 2 = 1; 0)3 = 2; 0)4 = 10, 0)5 = 15, Cjq = 20. 
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Fig. 2. Solutions of equation fl2^ with Cbp = 331, 35 and Ui = {0, 2; 1; 2; 10; 15; 20}, 
which correspond to Tcp = 27, 325; Tc^ = 30, 0255; 

1^3 = 32,1397; f}, = 41, 8727; = 46,1863; = 49, 9754. 


It is seen from £g.2 that the critical temperature grows with increasing phonon frequency ooq. 
The relations of critical temperatures Td/oJoi corresponding to the chosen parameter values are 
given in Table I. Table I suggests that the critical temperature of a Tl-bipolaron gas is always 
higher than that of ideal Bose-gas (IBG). It is also evident from fig.2 that an increase in the 
concentration of Tl-bipolarons n will lead to an increase in the critical temperature, while a 
gain in the electron mass m* to its decrease. For a) = 0 the results go over into the limit of 
IBG. In particular, fl22ll for a; = 0, yields the expression for the critical temperature of IBG: 

Te = 3, 3i;fV/V^e (23) 


It should be stressed, however, that fl2^ involves Mg = 2m*, rather than the bipolaron mass. 
This resolves the problem of the low temperature of condensation which arises both in the 
SRP theory and in the LRP theory in which expression fl25]l involves the bipolaron mass [26] - 
|29j . Another important result is that the critical temperature Tc for the parameter values 
considerably exceeds the gap energy coq. 

From (120]) . fl2TD follows that: 


N'{u) _ T3/2 fCj 
N ~~C^ '/'Vt 


(24) 


7Vo(c5) N\u) 

N N 

10 


(25) 













Fig.3 shows temperature dependencies of the number of supracondensate particles N' and the 
number of particles Nq occurring in the condensate for the above-indicated parameter values 



1-NVN 



Fig. 3. Temperature dependencies of the relative number of supracondensate particles N'/N 
and the particles occurring in the condensate Nq/N = 1 — N'/N for the parameter values Ui, 

given in Fig. 2. 

Fig.3 suggests that, as could be expected, the number of particles in the condensate grows as 
the gap Ui increases. 

The energy of a Tl-bipolaron gas E is determined by the expression; 

E = ^ rhnEn = EbpNo + ^ rhnEn (26) 

n=0,l,2,... 
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With the use of ffn|) . ffT^ . fl2B]) the specihc energy (i.e. the energy per one Tl-bipolaron) 
E{T) = E/Nu*, Ehp = Ehp/uj* will be: 


E{f) 


Ebp + -^E‘i/2 

^bp 


(¥) 


U ^ ^5/2 (v^) 

^ ^ (^) 


-^5/2 (a) 


J_ x^/^dx 
Jo 6^“^“ - 1 ’ 


where /i is determined from the equation: 


(27) 


(^) = (28) 

_ ^ f 0, f < f,; 

'* I A(n. r > f. 

Relation of jl with the chemical potential of the system /i is given by the formula fi = {fi — 
Ebp)/oj*. From fl27|) -fl25 ]) also follow expressions for the free energy: E = —2i?/3 and entropy 
5 = -dE/dT. 

Fig.4 illustrates temperature dependencies /S.E = E — Ebp for the above-indicated parameter 
values Ui- Break points on the curves AEi{T) correspond to the values of critical temperatures 

Te.. 


AE 
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Fig. 4. Temperature dependencies AE{T) = E{T) — Ehp 
for the parameter values a), presented in Fig. 2, 3. 


The dependencies obtained enable us to find the heat capacity of a Tl-bipolaron gas: C^(T) = 
dE/dT. With the use of fl7r|) for T < Tc we express C^{T) as: 

r.~,2 


C,{T) = 


f3/2 


2Cbp LT2' 


T/2 


+ 6 


" 3/2 




(29) 


r. / N 2 /■“ 1 dx 

vSe'+”-l 

Expression (12^ yields a well-known exponential dependence of the heat capacity at low tem¬ 
peratures ^ exp(—(Uo/T), caused by the availability of the energy gap coo- 

Fig.5 shows temperature dependencies of the heat capacity C^{T) for the above-indicated 
parameter values cDj. Table I lists the values of the heat capacity jumps: 


df 


dC,{f) 

df 


r=rc+o 


dC,{f) 

df 


f=Tc-0 


(30) 


at the transition points for the parameter values a)*. 

The dependencies obtained will enable one to find the latent heat of transition q = TS, 
where S is the entropy of supracondensate particles. At the point of transition this value is: 
q = 2T(.C^{Tc — 0)/3, where C'.y(T) is determined by formula fl2^ . For the above-indicated 
parameter values Wi, it is given in Table I. 
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Fig. 5. Temperature dependencies of the heat capacity for various 


values of the parameters cup uq = 

0 

0 

, = 25,2445; C 

viTcll) -- 

= 1,925 

UJi 

= 0,2 

; Tci ■ 

= 27,325; 

C,{fci 

-0) = 

2,162; 

Cvifcl 

+ 0) = 

1 ,868; 

UI 2 

= 1; 

P2 = 

30,0255; 

a{fc2 

-0) = 

2,465; 

Cv{fc2 

+ 0) = 

1,812; 

U 3 

= 2; 

fc3 = 

32,1397; 

Cy{fcs 

-0) = 

2,699; 

Cvits 

+ 0) = 

1,777; 

OJi - 

= 10; 

Tci = 

41,8727; 

CviTci • 

-0) = 

3,740; 

Cv{Tc4 

+ 0) = 

1,6779; 

(Us 

= 15; 

Tc5 = 

: 46,1863; 

Cv{Tc5 

-0) = 

4,181; 

Cv {Tc5 

+ 0) = 

1,651; 

ue 

= 20; 

Tc& = 

: 49, 9754; 

Cv (Tc6 

-0) = 

4,560; 

Cv{Tc% 

+ 0) = 

1,633. 


5 Current states of a Tl-bipolaron gas. 

In the foregoing we have considered equilibrium properties of a Tl-bipolaron gas. The for¬ 
mation of Bose-condensate per se does not mean that it has superconducting properties. To 
demonstrate such a possibility let us consider the total momentum of a Tl-bipolaron: 

V = Pi + P 2 + ka'^ttk, (31) 

where Pi and P 2 are the momenta of the hrst and second electron, respectively. It is easy to 
check that P commutes with Hamiltonian ([T]) and therefore is a constant value, i.e. c-number. 
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For this reason, to consider noneqnilibrium properties and, particnlarly cnrrent states we 
can nse generalization of Heizenberg transformation (jl]) snch that: 


Si{V) = exp 



ka 


k 


R. 


(32) 


The general expression for the fnnctional of the total energy of a Tl-bipolaron for P 7 ^ 0 is 
given in |3T] . Tl-bipolarons occnrring in condensed state have a common wave function for the 
whole condensate and do not thermalize in the state when P 7 ^ 0. In the supracondensate part, 
phonons whose wave vectors contribute into P fl3T|l are thermalized and their total momentum 
will be equal to zero. 

Hence, all the changes arising in considering the case of P 7 ^ 0 concern only the expression 
for the ground state energy which become dependent of P, while the spectrum of excited states 
flTT]) . flT^ remains unchanged. It follows that in passing over the critical point, the currentless 
state suddenly becomes current which is in agreement with the experiment. 


6 Comparison with the experiment. 

Fig. 4 shows typical dependencies of E{T). They suggest that at the point of transition the en¬ 
ergy is a continuous function of T. This means that the transition per se occurs without energy 
expenditure being a phase transition of the 2 -kind in complete agreement with the experiment. 
At the same time transition of Bose particles from a condensate state to a supracondensate 
one occurs with consumption of energy which is determined by the value q (§4, Table 1), de¬ 
termining the latent heat of transition of a Bose gas which makes it a phase transition of the 
1 -st kind. 

By way of example let us consider HTSC YBa 2 Cu^Oj with the temperature of transition 
90 -h 93K, volume of the unit cell 0,1734 • 10“^^ cm“^, concentration of holes n ~ 10^^ cm“^. 
According to estimates |13], Fermi energy is equal to: ep = 0,37 eV. Concentration of TI- 
bipolarons in YBa 2 Cu 307 is found from equation fl 2 ^ : 


= Ut) (33) 

n 

with Tc = 1, 6 . Table I lists the values of nbp^i for the values of parameters Ui given in §4. 
It follows from Table I that n^p^i « n. Hence, only a small part of charge carriers is in a 
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bipolaron state which justihes the approximation of a low-density Tl-bipolaron gas used by us. 
The energy levels of such Tl-bipolarons lie near Fermi surface and are described by the wave 
function: 

^(f) = e*^^>(r), (34) 

which leads to replacement of T, involved in (jH]), by: 

f = -2((p|A,|(p) + 24, (35) 

that is to reckoning of the energy from Fermi level (the last term in the right-hand side of fl35|) 
in dimensional units is equal to 26^?, where ep = h^k%/2m*). 

According to our approach, superconductivity arises when coupled states are formed. The 
condition for the formation of such states has the form: 

\Ebp\ > 2 \Ep\ , (36) 

where Ep is the energy of a Tl-polaron |3B]. Condition fl5B]) determines the value of a pseudogap: 

Ai = \Ehp — 2Ep\ (37) 

For Ai >> ojq the value of a pseudogap can greatly exceed both Tc, and the energy of the 
gap (i.e. cuo)- The expression for the spectrum E^p and Ep ffTT)) - ffT5|l suggests that the angular 
dependence Ai(fc) is completely determined by the symmetry of the isoenergetic surface of the 
phonon wave vector k. Earlier this conclusion was made by Bennet m, who proved that 
the main source of anisotropy of superconducting properties is the angular dependence of the 
phonon spectrum, though some contribution is also made by the anisotropy of Fermi surface. 

It follows from what has been said that formation of a pseudogap is a phase transition 
preceding the phase transition to the superconducting state. Recent experiments 051 also 
testify in favor of this statement. 

In paper [3l] correlation length for Tl-bipolarons was calculated. According to in 

HTSC materials its value can vary from several angstroms to several tens of angstroms, which 
is also in agreement with the experiment. 

Of special interest is to determine the characteristic energy of phonons responsible for the 
formation of Tl-bipolarons and superconducting properties of oxide ceramics. To do this let us 
compare the calculated values of the heat capacity jumps with experimental data. 
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As is known, in BCS theory a jnmp in the heat capacity is equal to: 


a - a 


a 


1,43 


where Cg is the heat capacity in the superconducting phase, and Cn is the heat capacity in the 
normal phase, and is independent of the parameters of the model Hamiltonian. As it follows 
from numerical calculations shown in £g.5 and in Table I, as distinct from the BCS theory, the 
value of the jump depends on the phonon frequency. Hence, the approach presented predicts 
the existence of an isotopy effect for the heat capacity jump. 



Fig.6. Comparison of the theoretical (solid line) and experimental |16] (broken line) 
dependencies in the region of the heat capacity jump. 

As it is seen from Fig.6, the heat capacity jump calculated theoretically (§4) coincides with the 
experimental value in YBa2Cu2,0^ [IH], for a; = 1,5, i.e. for cn = 7, 5 meV. This corresponds 
to the concentration of Tl-bipolarons equal to Ubp = 2, 6 • 10^® cm“^. Hence, in contrast to 
the widespread notion that in oxide ceramics superconductivity is determined by high-energy 
phonons (with energy 70 SOmeV ^\) actually, the superconductivity in HTSC materials 
should be determined by soft phonon modes. 

Notice that in calculations of the temperature of transition it was believed that the effective 
mass Me in equation fl22p is independent of the direction of the wave vector, i.e. isotropic case 
was dealt with. 

In the anisotropic case, choosing principal axes of vector fc, as coordinate axes, we will get 
the quantity {MexMeyM^zY^^ instead of the effective mass M^. In complex HTSC materials 
the values of effective masses lying in the plane of layers Mex, M^y are close in value. Assuming 
in this case Mg = Mgx = Mey = My, Mgz = Mx, we will get instead of Chp, determined by 
the value C^p = Cbpj^-, 7^ = Mx/My is the parameter of anisotropy. Hence anisotropy 
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of effective masses gives for the concentration the valne hf,p = 7 nfep. Therefore taking 
acconnt of anisotropy can an order of magnitnde and greater enhance the estimate of the 
concentration of Tl-bipolarons. If for YBa 2 CuzOi we take the estimate 7 ^ = 30 03. then for 
the concentration of Tl-bipolarons we will get: h^p = 1,4 • 10^® cm“^, which holds valid the 
general conclnsion: in the case nnder consideration only a small nnmber of charge carriers are 
in Tl-bipolaron state. The sitnation can change if the anisotropy parameter is very large. Thns, 
for example, in layered HTSC Bi-Sr-Ca-Cu-0 the anisotropy parameter is 7 > 100, accordingly, 
the concentration of Tl-bipolarons in these componnds can have the same order of magnitnde 
as the total concentration of charge carriers. 

Another important conclusion emerging from taking account of the anisotropy of effective 
masses is that the temperature of the transition depends not on Ubp and My individually, 
but on their relation which straightforwardly follows from ([ 22 ]). 

7 Essential generalizations of the theory. 

In the foregoing we considered the case of an ideal Tl-bipolaron gas. At small concentration of 
Tl-bipolarons their Coulomb interaction will be greatly screened which justihes the use of the 
model of an ideal gas. 

If the concentration of Tl-bipolarons is large (for example n = 10^^cm“^, as was believed 
in §4), then such Bose gas can no longer be considered to be ideal. Taking account of Coulomb 
interaction between bipolarons becomes necessary. For this purpose we can use Bogolyubov 
theory |l 8 ] for weakly imperfect Bose gas, which implies that the spectrum of elementary 
excitations with the momentum k will be determined by the expression: 

E{k) = \/k‘^u^{k) -|- {k‘^/2mBY (38) 

where u(k) = ^JnBV(k)/mB■, V{k) is the Fourier component of the potential of pairwise inter¬ 
action between charged bosons: V{k) = dvre^/eoA:^, and are the charge and mass of a 
boson, ns is the concentration of bosons. 

Whence it follows that: 

E{k) = {hwpY + {kP‘/2mBY, = y/dvre^ns/eomB (39) 

Up is the plasma frequency. According to (39), the spectrum of excitations of quasiparticles of 
an imperfect gas is characterized by a hnite energy gap which in the long-wavelength limit is 
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equal to Up. The same result can be arrived at if the Coulomb interaction between the electrons 
is considered as a result of electron-plasmon interaction. According to [IH], Hamiltonian of 
electron-plasmon interaction coincides in structure with Froehlich Hamiltonian which involves 
plasmon frequency instead of phonon frequency. This straightforwardly leads to the energy gap 
equal to the plasmon frequency cUp, if in (39) we put: = 2e, = 2m*, hb = n/2, where 

n is the electron concentration. Hence, for ojp < ujq the energy gap will be determined by Wp, 
while for Up > ujo h will still be determined by cuo- 

Actually, in real HTSC, there are not only phonon and plasmon branches, but also some 
other elementary excitations which can take part in electron pairing. An example is spin 
fluctuations. Generalization of the theory to the case of interaction with various brunches 
of excitations which will contribute into the ground state energy, the value of the gap and 
the dispersion law of a Tl-bipolaron is a topical problem. Obviously, taking account of this 
interaction will lead to an increase in the coupling energy of both Tl-bipolarons and Tl-polarons. 
Therefore, a priori, without any particular calculations one cannot say anything of how the 
condition of stability of Tl-bipolaron states determined by fl36|l will change. 

Another important problem is generalization of the theory to the case of intermediate cou¬ 
pling of electron-phonon interaction. Formally, the expression for the functional of the ground 
state energy of a Tl-bipolaron flT^ is valid for any value of the electron-phonon coupling con¬ 
stant. For this reason such a calculation will not change the spectrum of a Tl-bipolaron, 
however it will alter the criteria of fulhllment of the conditions of a Tl-bipolaron stability (I5B]) . 

8 Conclusive remarks. 

In this paper we have presented conclusions emerging from consistent translation-invariant 
consideration of EPI. It implies that, pairing of electrons, for any coupling constant, leads 
to a concept of Tl-polarons and Tl-bipolarons. Being bosons, Tl-bipolarons can experience 
Bose condensation leading to superconductivity. Let us list the main results following from 
this approach. First and foremost the theory resolves the problem of the great value of the 
bipolaron effective mass (§4). As a consequence, formal limitations on the value of the critical 
temperature of the transition are eliminated too. The theory quantitatively explains such 
thermodynamic properties of HTSC-conductors as availability (§4) and value (§6) of the jump 
in the heat capacity lacking in the theory of Bose condensation of an ideal gas. The theory also 
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gives an insight into the occnrrence of a great ratio between the width of the psendogap and 
(§6). It acconnts for the small valne of the correlation length [M] and explains the availability 
of a gap (§3) and a psendogap (§6) in HTSC materials. The angular dependence of the gap 
and psendogap gets a natural explanation (§6). 

Accordingly, isotopic effect automatically follows from expression 022 p . where the phonon 
frequency ujq acts as a gap. The conclusion of the dependence of the temperature of the 
transition on the relation nfep/My (§6) correlates with Uemura law universal for all HTSC 
materials which implies that the temperature of the transition relates to the concentration of 
charge carriers divided by its effective mass [50] . 

Application of the theory to ID and 2D systems leads to qualitatively new results since the 
occurrence of a gap in the Tl-bipolaron spectrum automatically removes divergences at small 
momenta, inherent in the theory of ideal Bose gas. These problems will be considered by the 
author in a special paper. The case of ID-polaron was discussed in [51] . 

The work was done with the support from the RFBR, project N 13-07-00256. 


9 Appendix. Remarks on the notation. 

Function F^i 2 {cx) is called a polylogarithm = Li 3 / 2 (e“"), in mathematics this is the function 
PolyLog, therefore the function in fl2^ will be: = f^/‘^PolyLog 3/2, . 

In the general case, function PolyLog of the order of s is determined as: 

1 7 “ P- ' 

i^-1' 


PolyLog [s, e "] = 


-dt 


where F(s) is a gamma function: F(l/2) = F(3/2) = \/^/2, F(5/2) = 3i/F/4. 

Accordingly, the functions R1/2, T3/2, F5/2 occurring in the text will be: 

Fi /2 = 2PolyLog 1/2, 

F 3/2 = PolyLog 3/2, 

-P 5/2 = ^/2PolyLog 5/2, 

References 


[ 1 ] L.P. Pitaevskii, E.M. Lifshitz, Statistical Physics, Part 2: Theory of the Condensed State 
(1st ed.). Butterworth-Heinemann. ISBN 978-0-7506-2636-1. (1980). 


21 












[2] J.G.Bednorz, K.A.Miiller, Z.Phys.B 64, 189 (1986). 

[3] Ginzburg V.L. Phys. Usp. 43, 573 (2000); DOB 10.1070/PU2000v043n06ABEH000779 . 

[4] S.L.Kakani, S.Kakani, Superconductivity, Anshan, Kent, UK (2009). 

[5] T.Tohyama, Jap. J.Appl.Phys. 51, Article ID 010004 (2012). 

[6] S.Kruchinin, H.Nagao, S.Aono, Modern Aspects of superconductivity. Theory of supercon¬ 
ductivity, World Sci., River Edge, NJ USA (2011). 

[7] K.P.Sinha, S.L.Kakani, PNAS 72, 153 (2002). 

[8] K.H.Benneman, J.B.Ketterson, Superconductivity: Gonventional and Unconventional Su¬ 
perconductors 1-2, Springer, New-York, NY, UK (2008). 

[9] J.R.Schrieffer, Theory of Superconductivity Westview Press, Oxford, UK (1999). 

[10] A.S.Alexandrov, Theory of Superconductivity from weak to strong oupling, lOP, publish¬ 
ing, Bristol, UK (2003). 

[11] N.M.Plakida, High Temperature Guprate Superconductors: Experiment, Theory and Ap¬ 
plications, Springer, Heidelberg, Germany (2010). 

[12] I.Askerzade, Unconventional Superconductors 153, of Springer Series in Material Science, 
Springer, Berlin, Germany, (2012). 

[13] O.Gunnarsson, O.Rosch, Journ. of Phys. 20, Article ID 043201 (2008). 

[14] T.Moria, K.Ueda, Adv. Phys. 49, 556 (2000). 

[15] D.Manske, Theory of Unconventional Superconductors, Springer, Heidelberg, Germany 
(2004). 

[16] J.Bardeen, L.N.Gooper, J.R.Schrieffer, Phys. Rev. 108, 1175 (1957). 

[17] A.S.Alexandrov, N.Mott, Polarons & Bipolarons, World Sci. Pub. GO Inc, Singapore 
(1996). 

[18] Alexandrov A.S., Krebs A.B. Sov. Phys. Usp. 35, 345 (1992). 


22 



[19] P.W.Anderson, Physics Word 8, 37 (1995). 

[20] N.F.Mott, Physics World 9, 16 (1996). 

[21] P.W.Anderson, Physics World 9, 16 (1996). 

[22] E.V.L.Mello, J.Ranninger, Phys. Rev. B 55, 14872 (1997). 

[23] Yu.A.Firsov, V.V.Kabanov, E.K.Kudinov, A.S.Alexandrov, Phys. Rev. 59, 12132 (1999). 

[24] E.V.L.Mello, J.Ranninger, Phys. Rev. B 59, 12135 (1999). 

[25] R.A.Ogg, Phys. Rev. 70, 93 (1946). 

[26] Vinetskii V.L., Pashitskii E.A., Ukrainian Journal of Physics 20, 338 (1975). 

[27] Pashitskii E.A., Vinetskii V.L. JETP Lett. 46 (SI) S104 (1987). 

[28] D.Emin Phys. Rev. Lett. 62, 1544 (1989). 

[29] Vinetskii V.L., Kashirina N.L, Pashitskii E.A., Ukrainian Journal of Physics 37, 76 (1992). 

[30] M.R.Schafroth, Phys. Rev. 96, 1149 (1954). 

[31] Tulub A.V. Vestnik Leningrad Univ. 22, 104 (1960). 

[32] Tulub A.V. Sov. Phys. JETP 14, 1301 (1962). http://www.jetp.ac.rU/cgi-bin/e/index/r/41/6/pl828? 

[33] Lakhno V.D. JETP 110, 811 (2010). 

[34] V.D.Lakhno, Sol. St. Comm. 152, 621 (2012). 

[35] Kashirina N.L, Lakhno V.D., Tulub A.V. JETP 114, 867 (2012). 

[36] Lakhno V.D. JETP 116, 892 (2013). 

[37] Lakhno V.D. Phys. Usp. 58, 295 (2015). 

[38] W.Heisenberg, Zs. fiir Physik 65, 4 (1930). 

[39] T.Lee, F.Low, D.Pines, Phys. Rev. 90, 297 (1953). 

[40] Tyablikov S.V. Methods of Quantum Theory of Magnetizm (M.: Nauka, 1975) 


23 


[41] S.J. Miyake, in Polarons and Applications, V.D. Lakhno (Ed.), Wiley, Chichester, (1994), 
514 p. 

[42] Levenson O.B. Rashba E.I., Usp. Fiz. Nauk 111, 683, (1973); Sov. Phys. Usp. 16, 892, 
(1974) 

[43] Gorkov L.P., Kopnin N.B. Sov. Phys. Usp. 31, 850 (1988). DOI; 

10.1070/PU1988v031n09ABEH005623 

[44] A.J.Bennet, Phys. Rev. 140, A1902 (1966). 

[45] Rui-Hua He, M.Hashimoto et ah. Science 331, 1579 (2011). 

[46] N.Overend, M.A.Howson, I.D.Lawrie, Phys. Rev. Lett. 72, 3238 (1994). 

[47] A.Marouchkine, Room-Temperature Superconductivity, Cambridge Int. Sci. Pubh, Cam¬ 
bridge (2004). 

[48] Bogolyubov N.N. J. Phys. USSR 11, 23 (1947). 

[49] V.D.Lakhno, in Polarons and Applications ed. by V.D.Lakhno, Wiley, Chichester (1994), 
p.514. 

[50] Y.J.Uemura et ah, Phys. Rev. Lett. 62, 2317 (1989); 66, 2665 (1991). 

[51] V.D.Lakhno, Progr. Theor. Exp. Phys. 073101 (2014), doi: 10.1093/ptep/ptu075 


24 



Table I. Calculated characteristics of Bose-gas of Tl-bipolarons with concentration n = 10^^ 

cm”^. 


i 

0 

1 

2 

3 

4 

5 

6 

Ui 

0 

0,2 

1 

2 

10 

15 

20 

Tci /^oi 

oo 

136,6 

30,03 

16,07 

4,187 

3,079 

2,499 

qi/Tci 

1,28 

1,44 

1,64 

1,80 

2,49 

2,79 

3,04 

-^{dC^JdT) 

0,114 

0,119 

0,124 

0,128 

0,143 

0,148 

0,152 

Cv,i{Tc — 0 ) 

1,925 

2,162 

2,465 

2,699 

3,740 

4,181 

4,560 

(C, - Cn)/Cn 

0 

0,1574 

0,3604 

0,519 

1,229 

1,532 

1,792 

rihpi-cvsY 

le-io^^ 

9,4-10i« 

o 

00 

2 , 0 - 101 ^ 

iggiH 

5,2-10i4 

2,3-101^ 


oji = oji/oj*, u* = SmeV, ca* is the energy of an optical phonon; Tc^ is the critical temperature 

of the transition, g* is the latent heat of the transition from condensate to supracondensate state; 

—AidCvi/dT) = dCvi/dT _ _ — dCvi/dT _ is the jump in heat capacity during SC 

’ ’ f=f^^+o ’ T='h.-o 

transition, T = T/u* ; Cy^iiTc — 0) is the heat capacity in the SC phase at the critical point; 
Cs = C^{Tc — 0), = C^{Tc + 0); The calculations are carries out for the concentration of 

Tl-bipolarons n = 10^^ cm“^ and the effective mass of a band electron m* = rriQ. 

The table also lists the values of concentrations of Tl-bipolarons rihpi for HTSC YBa2CuzOT, 
based on the experimental value of the transition temperature = 93i^ (§ 6 ). 
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